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Abstract. Driven lattice gases as the ASEP are useful tools for the modeling of 
various stochastic transport processes carried out by self-driven particles, such as 
^ ' molecular motors or vehicles in road traffic. Often these processes take place in one- 

\l , dimensional systems offering several tracks to the particles, and in many cases the 

particles are able to change track with a given rate. In this work we consider the case 
of strong coupling where the hopping rate along the tracks and the exchange rates are 
of the same order, and show how a phenomenological approach based on a domain wall 
r^ ' theory can describe the dynamics of the system. In particular, the domain walls on 

^3 . the different tracks form pairs, whose dynamics dominate the behavior of the system. 
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1. Introduction 

The dynamics of particles on microscopic scales is generically non-deterministic since the 
motion of the particles is influenced by the surrounding medium. The stochastic motion 
can either be passive as for Brownian particles or active, i.e. particles are self-driven by 
processes consuming energy. 

In biological systems molecular motors are example of active particles, which 
combine directed and stochastic motion [1]. The molecular motors are proteins carrying 
out a large variety of intracellular transport processes. Their active dynamics is 
characterised by a directed stepwise motion along the cytoskeleton, a network of one- 
dimensional filaments. The energy of this process is provided by the hydrolysis of ATP. 
A detailed understanding of the motor dynamics is of obvious biological interest. 

Next to the microscopic systems self-driven particles exist also on larger scales, e.g. 
models describing vehicular traffic |2l |3]. The stochastic nature of the motion is in this 
case not related to thermal fluctuations, but rather to the integration of uncontrolled 
degrees of freedom - the distraction of a driver, its psychological state, the mechanical 
reactions of the car, etc. In the same way, pedestrian or gregarious animal traffic can be 
classified in the same family of interacting self-driven particles described by stochastic 
processes (see e.g. [1]). 

From a theoretical point of view, it is of interest to extract the generic non- 
equilibrium behaviour of these systems. Systems of interacting self-driven particles 
show particularly interesting behavior if open boundary conditions are applied. Then 
one observes so-called boundary induced phase transitions, i.e. the boundary particle- 
reservoirs determine the bulk-properties of the system [5j . The transport properties and 
particle distributions are the result of a subtle interplay between boundary reservoirs and 
bulk properties of the system, which has been analysed in great detail. The paradigmatic 
model for self-driven particles is the so-called asymmetric exclusion process (ASEP) 
[6l [7]. In the case of the ASEP particles perform a directed random- walk on a one- 
dimensional lattice. Particles interact via steric exclusion, i.e. lattice sites are either 
empty or occupied by a single particle. By using recursion relations the stationary 
properties of this model have been exactly calculated. 

Although an exact solution of the model exists for the stationary state, it is 
instructive to exploit the picture of competing boundary reservoirs |8j directly, as it 
can be extended to non- stationary regimes [9] : Each particle reservoir tries to impose a 
domain of constant density. Then the bulk properties of the system are determined by 
the dynamics of the shock, i.e. the position where one observes the sudden transition 
from one domain to the other. The dynamics of the shock can be described as a (biased) 
random walk with constant rates, where the left and right hopping rates depend on the 
flow realized in the two domains. It is important to notice that the hopping rates of 
the shock do not depend on its position, a property which is closely related to the 
conservation of particles in the bulk of the system [lOl [HI [121 (El [H] ■ 

Variants of the ASEP have been used in order to model a number of important 
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real empirical systems, including vehicular and intracellular traffic. In order to apply 
particle hopping models to vehicular traffic one has to introduce at least different particle 
velocities and a finite reaction time of the drivers via a parallel update scheme [T5] . 
This reaction time can be enhanced by making the hopping probabilities dependent 
on the velocity [161 [E]- While traffic models based on cellular automata give today 
a satisfying description of intra- lane dynamics [18], there is a growing interest in the 
traffic community to understand better the lateral interactions and the dynamics of 
lane-changes. 

In the present work we analyse a generalized ASEP which is defined on a two-lane 
system. We consider uni-directional motion, i.e. particles move in the same direction 
on both lanes. In the bulk, the particle number is conserved globally, but this is not 
the case for each separate lane, since lane-changes are allowed. This setup is relevant 
for vehicular traffic on multilane roads, pedestrian traffic and intracellular transport. 
A model of this kind has been studied recently by Reichenbach et al. [TU [201 EI] who 
concentrated on the weak coupling limit, which is relevant for processive molecular 
motors. In this limit the small amplitude of the exchange rates implies a continuous 
flow proffie for each track. 

Here we consider the system for arbitrarily strong lane-changing rates, i.e. exchange 
rates that are at least of the same order as the hop rates, which is relevant in case of 
pedestrian and vehicular traffic, where frequent lane-changes are observed. In this case 
it is possible to observe not only a discontinuous density profile but also a discontinuous 
flow profile. We show here that the system can still be characterized by a domain wall 
approach. In particular we concentrate on the dynamics and synchronisation of the 
shocks in the system, and show how they dominate the system behavior. 

The paper is organised as follows. First a review of the results obtained in the 
case of weak coupling will be presented. As preliminaries to the strong coupling case, 
a notion of adaptation between the lanes will be defined, and the possibility to extend 
the use of second class particles (SCP) [221 123] to two-lane systems will be discussed. 

The dynamics of shocks under strong coupling will be exemplified on a special case, 
illustrating how shocks on both lanes are coupled and how their joint motion can be 
predicted in a domain wall picture. We also discuss the inner structure of the shock 
pair and the relevance of our results for the phenomena taking place in strongly coupled 
driven systems in general. 

2. A tvi^o-Iane model 

2.1. Definition of the model 

In this work we consider a L x 2 lattice with open boundary conditions. Each lattice 
site may either be empty or occupied with a single particle. In the bulk particles hop 
from site i — )■ i -|- 1, where 1 < i < (-^ — 1), with rate p^{p^) in upper (lower) lane if site 
i + 1 empty. In the remaining of the paper, we shall always use p^ = p^ = 1. Particles 
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Figure 1. Schematic illustration of the particle dynamics. 



may also change from the upper to the lower lane with rate uj^ and with rate w^ in the 
inverse direction. 

The asymmetry between the vertical couplings is defined by u^ = Acj^. The density 
of particles is controlled via left and right particle reservoirs which are coupled to the 
chains. Particles enter the system at site 1 with rate a^ (upper lane) or rate a^ (lower 
lane), if the first site is empty. The outfiow from the chain is controlled by the rates /3^ 
(upper lane) or (3^ (lower lane). 

We shall denote by j^ and j^ the horizontal fiuxes of particles respectively along 
the upper and lower lane, while h^ and h-^ stand for the vertical fiuxes from the lower 
to the upper (resp. upper to lower) lanes and the net vertical fiux is /i"^* = h^ — h^. 



2.2. Mean-field equations 

The model introduced in the previous section is similar to the one in [20], and we recall 
how the mean-field equations can be derived. 

The mean-field expressions for these fiuxes are 



fix + 2^ 
1, 



p'[x^ 



[I - p^ix + 1 



j\x + -)= p^{x) [1 - p^{x + 1)] 



h^x) 
h\x) 



u^p'^{x) [l — p^{x)] 
tu^p'^ix) [l - p^{x)] 



In a stationary state, mass conservation implies that 

fix + 2'' ~ -^^^^ ~ 2'' " ^^^^"^ " ^^^^"^ 
fix + 2^ ~ ^^^^ ~ 2'' " ^^*'^'' ~ ^^*-^'* 



(1) 

(2) 

(3) 
(4) 

(5) 
(6) 



We now consider the limit L ^ 1 and introduce a new space variable y = x/L. 
Then the mean-field equations read 



(1 - 2p^) dyp'^ = dyf = u^Lp^ (1 - p^) - uj^Lp^ (1 - P^) 
(1 - 2p^) dyp^ = dyj^ = u^Lp^ (1 - p^) - co^Lp^ (l - p^) 



(7) 
(8) 
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i.e. we obtained coupled differential equations for the densities. Obviously two 
independent ASEPs are obtained if the coupling rates u^ and w^ vanish. 

These equations will be valid only in the slowly varying regions of the density 
profile. 

3. Overview of useful notions 

3.1. Domain wall picture for one-lane ASEP 

The stationary and the transient state of driven lattice gases with open boundary 
conditions are characterized by a competition between the capacity of the boundary 
reservoirs and the capacity of the chain. This feature of driven lattice gases has been 
used in order to develop an effective theory for the macroscopic properties of open 
systems [8]. This section can be ignored by the reader already familiar with these 
notions. 

The basic idea of the domain-wall-picture is most simply understood for one- 
dimensional driven lattice gases with mass conservation and a single particle-species. 
The most prominent example of this type of models is the TASEP (Totally asymmetric 
exclusion process) which we will consider for concreteness. In this case the left reservoir 
of capacity a < a^ax m establishes a low density domain of density pi (a) at the left of 
the system. The current in the low density domain is given by j{a) (if j{a) < jmax) 
where jmax is the capacity of the chain. The low density domain competes with a high 
density domain of density Pr-(/3) which is controlled by the right reservoir of capacity /3 
if (j(/3) < jmax)- Both domains are separated by a so-called domain wall or shock which 
is located at Xs{t) \§. In general we have pi{a) ^ Pr{(^) such that the density profile is 
discontinuous at the shock position Xgit). 

In order to characterize the properties of the system it is crucial to describe the 
dynamics of the domain wall. The dynamics of the domain wall can be illustrated by 
considering a half-open chain first, i.e. an ASEP with (/3 = 0). We also assume that 
the chain is initially empty, and thus that the domain wall is initially positioned at the 
right boundary. After introducing A^ particles to the system, the domain wall is shifted 
A^-sites to the left assuming that a ^ p. As A^ = j(a)At the average velocity of the 
domain wall is simply given by j{a). 

For arbitrary values of a, /3 one has to consider the displacement for general values 
of the densities inside the coexisting domains. In this case A^ particles correspond to a 
displacement Ax = N/{pr{f3) — pi{oi)) and the average velocity of the domain wall is 
given by Vr = j{a)/{pr{f3) — pi{a)). Analogous the motion of the wall towards the exit 
has the velocity vi = j {(3) / {pr{(3) — pi{a)). 

Due to the stochastic nature of the coupling between boundary and chain the motion 
of the domain wall is well described as a random walk with right and left hopping rates 

I The boundary rates are given as multiples of p. a„iax = 1/2 for the ASEP. 

§ Note that the discussion refers to the continuum limit of the model. In discrete space we can only 

distinguish between empty and occupied sites. 
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Vr and Vi- From this picture, it is possible to derive (with a very good agreement 
with direct simulations) the density profile, but also more subtle quantities such as the 
fluctuations of the number of particles in the system [9] or the largest relaxation time 
in the system [SI [2l] . 

This picture is valid for general one dimensional systems with particle conservation 
in the bulk of the system. For systems that are coupled to bulk reservoirs, i.e. for 
which the particle number is not conserved even in the bulk, the above picture has 
to be modified. Due to particle exchange the flow inside the domains is no longer a 
conserved quantity. For such systems the hopping rates of the domain wall are position 
dependent. This implies for example that the shock can be localized in the bulk at the 
position < Xs < L if the bias of the walker changes its sign at Xg [HI [121 US]- The 
present model combines global particle conservation with particle exchange between the 
two lanes. Therefore one expects to observe localization effects as well as the coupling 
between the motion of the two domain walls, one on each lane, in the system. 

3.2. The case of weak coupling - a review 

Different variants of the present model have been considered recently where the coupling 
of the lanes has been realized in different ways. The first class of models is characterized 
by hopping rules which respect mass conservation in each lane in the bulk [25l |26l ET] . 
In these models the coupling between the two lanes has been realized via hopping rates 
depending on the configuration on the other lane, but there are no particle exchanges 
between the lanes. Models with non-vanishing exchange rates have also been considered 
before [28]. Pronina and Kolomeisky [29] studied the stationary state of a two-lane 
system for symmetric and asymmetric lane changing rules which are coupled to particle 
reservoirs. The in- and output rates are the same for the two chains. As a consequence 
shocks are localized at the same position for the two lanes and the structure of the phase 
diagram can be obtained by a mean-field analysis. 

More general capacities of the boundary reservoirs have been considered in a series 
of papers by Reichenbach et al. [201 EI]- In their work the present model has also been 
discussed in terms of a spin representation. Their focus was on establishing the phase 
diagram and analysing the systems behaviour for weakly coupled chains, i.e. the case 
where the exchange rates scale as 1/ L. From a technical point of view this assumption 
implies that that the current between the two lanes is also of the order 1/L. Therefore 
the flow profiles are continuous at arbitrary positions in the system. This is the main 
difference to our work since here both the exchange rates as well as the density difference 
between upper and lower lane are of C(l). 

4. Adaptation between lanes 

Focussing on the domain wall dynamics in the strong coupling limit we choose the model 
parameters such that no additional boundary effects inside coexisting domains appear. 
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as this will be explained in this section. 

For a < 1/2 and (5 < 1/2, as we have seen in the previous section, a one-lane 
system can be described by the domain wall picture: a domain wall (or density jump) 
separates two domains of constant density. Now the question is: What happens when 
two flat density profiles on two different lanes are coupled (a; 7^ 0)? In the general case, 
a vertical flux immediately sets in. Assuming a factorized structure on each lane, the 
vertical fluxes just after the onset of the coupling are 

h^x) = u;V(x) [1 - p^{x)] (9) 

h^{x) = a;V(a;) [l - p\x)] . (10) 

It may happen that p^ and p^ are such that the vertical net flux vanishes. In this 
case the densities are conserved on each lane, and the flat proflles are maintained. We 
deflne density relations which conserve the densities in both lanes as adapted. With a 
relation u^ = \u^, the condition for the upper and lower densities to be adapted reads 

by deflnition of fx. We also notice that 

f\/,{X) = f-\X). (12) 

In the general case, i.e. if the initial densities are not adapted, a non zero vertical 
flux occurs for non-vanishing exchange rates. In the bulk of the system a new local 
equilibrium between the lanes is reached quite rapidly (typically a few time steps times 
the inverse exchange rates for the parameters of table [I]). The stationary bulk densities 
are adapted and the vertical net flux vanishes. 

At the boundaries, each reservoir tries to impose a density in the system. As the 
upper and lower lane reservoirs are decoupled, the particle densities near the boundaries 
are not automatically adapted and it is possible to observe a stationary net flow between 
the two lanes in the vicinity of the entrance or exit. If the entrance rate on one given 
lane is a^ (or a'^), it is equivalent to say that the entrance reservoir tries to impose a 
density a^ (resp. a^) in the system. If the densities a^ and a^ of the entrance reservoirs 
fulflU the relation (fTTj) . then we say that the entrance rates are adapted. The condition 
reads 

°'= l + (A-l)„t -M°^) («) 

In the same way, the exit reservoirs try to impose the densities 1-/3^ and 1 — P^ 
respectively on the upper and lower lane. Therefore the condition for the exit rates to 
be adapted is 

1-3^ = Mijl^!) ^ ^ (1 _ 0t) (14) 

^ l + (A-l)(l-/3t) ^^^ '' ' ^ ^ 



or equivalently 



P'- x-(x-,)p, -f^'^(^'^' (15) 
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parameter 


value 


parameter 


value 


at 


1/3 


/3^ 


0.176 


oc^ 


0.2 


/3^ 


0.3 


UJ^ 


0.2 


u^ 


Xu^ 


X 


1 

2 







Table 1. Set (I) of parameters. Boundary conditions are adapted on both sides. 

If boundary rates (entrance or exit) are not adapted a boundary layer forms where a 
non zero net flux between the lanes transfers from the reservoir densities to the adapted 
bulk densities (see an example in [30]). 

In the remaining of the paper, we choose only adapted boundary conditions, in 
order to avoid these boundary layers, and to focus on the role of domain walls. The case 
of non adapted boundary conditions will be discussed in section O Note that even if 
entrance and exit boundary rates are adapted, the adapted densities that the entrance 
reservoirs try to impose in the system are not necessarily the same as the adapted 
densities imposed by the exit reservoirs, hence one observes the formation of domain 
walls in between. 



5. Second class particles in two-lane systems 

In the case of one-lane ASEP systems with mass conservation, a useful tool for tracking 
shocks is the introduction of second class particles (SCP) [311 132]. A SCP is a passive 
tracer particle which does not modify the dynamics of the ordinary (or first class) 
particles. Its dynamics is defined in a way that, in presence of a discontinuity in the 
density profile of first class particles, i.e. a shock, the SCP will locate itself at the shock 
position (at least, for a certain class of shocks, as this will be explained below). Thus 
following the motion of the SCP allows to track the dynamics of the shock. 

The problem in using the SCP for this model is due to the fact that mass is not 
conserved for each lane. Therefore we have to modify the dynamics of the SCP in order 
to keep its ability to trace shocks also for two-lane systems. 



5.1. Definition of the dynamics for a second class particle in a two-lane TASEP 
The intra lane evolution rules are, as for one-lane systems: 

10 — ^01 
20 — ^02 
12 — ^21 

where 0, 1, and 2 stand for empty site, first class particle, and second class particle, 
respectively. The updates of the SCP are carried out with the hopping rates p^ and p^ 
of the first class particles. This rule implies that the SCP moves synchronously with 
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the shock [21] if the number of particles is conserved on each lane. For the model 
under consideration this is not the case as (first class) particles are exchanged with the 
boundary reservoirs and between the lanes. When a first class particle changes lane 
and the target site is occupied by a SCP one has to make sure that the SCP behaves 
like a hole. A simple exchange of the positions of first and second class particle is not 
possible since the SCP cannot change lane. Therefore, non-local moves of the SCP have 
to be introduced, as the neighboring sites of the SCP are not necessarily empty. In 
principle several rules are possible. In figure |2] we illustrate our choice of the non-local 
SCP dynamics. This rule ensures that the SCP does not change the lane and follows 
the domain wall position in the same lane. 

i-|7 I 6 I 5 I 3 I 1 |Q|2 |4 I , , , 



I I I I I I 1^1 I I I I I 

Figure 2. Schematic illustration of the second class particle (white disk) dynamics, 
when a first class particle (black disk) arrives from the other lane. Neighboring sites 
are explored in the order indicated by the figures, until an empty site is found. Then 
the second class particle moves to this empty site. 

When a first class particle changes lane and arrives on a site occupied by a SCP, 
the neighboring sites of the SCP are explored in the order that is indicated in figure |2j 
The SCP has to jump to the first empty site that is found. 

We have chosen an asymmetric move of the SCP in order to avoid the dynamical 
trapping of the SCP in the high density region to the right of the domain wall. Still, it 
has to be symmetric enough to avoid a too strong drift towards the left. 

5.2. Mean field calculation of the drift velocity of a SCP in a constant density 
background 

In this section, the drift velocity of a SCP placed in a flat density profile will be calculated 
(the term "density" always refers to a density of first class particles). We note p = 1 —p. 
In a mean-field approximation, the drift velocity of a SCP located on the upper lane 
would be 



4cp = (+i)p^ + (-i)p^ 



+ U^p^ 



-p^ + pV^ - 2p^'p^ + 2p^'p^ - 3p^V^ - V'p^ 



,t_,.,ta. 



2pi -w'p- 



1 - 3p^ + 5p^' - 7p^^ + 9p-^' - 4p^^ 



1-pt 

If the upper and lower densities are adapted, using (TTT]) . Vg(jp is a function only of p^ 
and u'^. In the special case of uncoupled lanes {u^^ = 0), the drift velocity v^^^p vanishes 
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exactly for p^ = 1/2. It is positive for a lower density and negative for a higher one. 
Thus the SCP locates itself on shocks that separate a left low density (<l/2) domain 
and a right high density (>l/2) domain. Note that shocks that do not go through the 
density 1/2 cannot be tracked. When the coupling is non zero, the density p^ for which 
Vg(jp vanishes is given by an implicit expression that can be solved numerically. The 
upper SCP will thus locate itself on shocks separating a left domain with density less 
than p^ and a right domain with density greater than p^ . 

In the same way, the density p^* for which Vg(jp vanishes is defined implicitely by 

[l - 3p^ + 5p^^ - 7p^^ + 9p^^ - 4p^^ 

= 1 - 2pJ- - Au;Vi/a(p^)^ , T (16) 

1 — p^ 

and thus the densities at which SCPs localize may be slightly different on the upper and 

lower lanes - though this difference turns out to be negligeable for the couplings that 

we are considering here. 

In particular, for the set of parameters of table [1] that will be used later, we find 

that p^ = 0.481 and p^ = 0.480. This is not too far from the one-lane value 0.5. Note 

however that the value of p^ depends on the coupling strength. Indeed, this is one 

reason why we have chosen the coupling constants of tables [1] and O For a too strong 

coupling, the value of p^ could be too close to the density of the right domain to ensure 

properly the tracking of the shock. 

6. The physics of strong coupling 

In order to explain the main phenomena involved in strongly coupled two-lane systems, 
we selected one representative set of parameters, given in[Tl The influence of the different 
model parameters is discussed in the text below. 

The influence of the shock is nicely illustrated in the transient regime. We therefore 
discuss behaviour of the model in terms of the different steps towards the stationary 
state. 

6.1. Initializing the system: Decoupled lanes 

We choose an initial condition such that the upper lane is in the high density (controlled 
by the exit rate) regime and the lower lane in the low density (controlled by the 
entrance rate) regime. More precisely for the entrance and exit parameters of table 
m and decoupled lanes, the bulk density in the upper lane (resp. lower lane) is given 
Plik = 1 - /3^ = 0.824 (resp. pL, = a^ = 0.2). 

The parameters are chosen such that the domain wall picture is valid, i.e. the 
entrance and exit rates are chosen such that they do not exceed the capacity of the 
chain. The domain wall is localized towards the entrance in the upper lane, and towards 
the exit in the lower lane. In order to track the dynamics of the domain walls we place 
one SCP on each lane. The decoupled system including the SCP is simulated until a 
stationary state is reached. 
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6.2. Early stage - local adaptation 

At t = 0, some coupling is added between the lanes (see parameters in table [1]). As 
outlined in section HJ this implies that a vertical flux immediately sets in. This is 
illustrated in Figure (|3] - a) showing the evolution of the density profile in the first 
20 time steps. In the bulk, particles are exchanged between the lanes until a local 
equilibrium is obtained, i.e. until the net flux between the lanes vanishes. This implies 
that at the end of this stage, the adaptation relation ( ITTl) is valid, i.e p^ = fx{pl). The 
actual values of the bulk densities can be obtained from local mass conservation: 

pl + pi = l-P^ + cy^ = Po (17) 

where the r.h.s. is the sum of bulk densities for uncoupled lanes and the l.h.s. is the 
sum of bulk densities 20 time steps after coupling is activated. 

From these two relations, we find that the bulk densities after 20 time steps should 
be given by 



t _ Po _ A + 1 
^* ~ 2 2(A - 1) 



A + 1 

Po 



-A^ (^* 



^A-1 

pi = Po-pI (19) 

where the '+' sign is valid if A > 1 and '-' if A < 1 (for the special case A = 1 one 
gets pi = pi = Po/2). In our case, A = 1/2, and the predicted bulk density values 
are pi = 0.598 and pi = 0.426, which are to be compared with the numerical values 
obtained in the simulations pi = 0.599 ± 0.002 and pi = 0.424 ± 0.002 (see Fig. [3]-a). 

6.3. Shock dynamics 

At the end of the local adaptation stage, a new domain wall appears on the right of 
the upper lane (resp. on the left of the lower lane ), to separate the new bulk domain 
from the exit driven (resp. entrance driven) domain. There are now four shocks in the 
system, located at both ends of each lane. 

The second class particles allow to track the shocks which are at t = on the left 
of the system in the upper lane and on the right of the system in the lower lane. The 
two other shocks cannot be traced by a SCP since they do not go through the density 
1/2. 

The two shocks on the right (resp. left) of the system are immediately coupled, 
and start to move in pairs towards the center of the system. Their drift velocity can be 
computed analytically (see next subsection). While the vertical flux is vanishing within 
each domain, it is nonzero between the two shocks of a given pair. This gives another 
signature of the motion of the pair, which can be compared with the SCP tracking. 

Eventually both shock pairs merge on each lane (see figures [3]-b and[3]-e). For our 
setup the two shock pairs merge in the bulk of the system around time 2300. The 
resulting pair of walls move towards the left boundary since the global capacity of the 
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Figure 3. Simulation results obtained for the parameters of tableland for a system 
size L — 1000. The left figures (a-c) show the density profiles at different times, while 
the right figures (d-f ) show the SCP locations in the corresponding time intervals. For 
all figures, solid (resp. dashed) lines refer to the upper (resp. lower) lane. Fig (a) : 
density profiles at time (most upper and lower curves, black online) and 20 (red 
online). Fig (b) : density profiles at time 20 (black, red online) and time 2304 (grey, 
green online). Fig (c) : density profiles at time 2304 (grey, green online) and time 
30000 (black, blue online). In each figure, light grey lines indicate 20 intermediate 
density profiles. All data have been averaged over 10^ trajectories. Error bars are 
given for the mean SCP trajectories; they are smaller than the line thickness and can 
hardly be seen. 



Shock- dynamics of two-lane driven lattice gases 13 

left reservoirs exceeds the capacity of the right reservoirs. In the final stationary state, 
it stays localized near the entrance (figures |3]-c and[3]-f). 

In figure [3] the time-dependent average position of the SCP is given. Actually, 
the SCP trajectories vary from one realization to another one, since it depends on the 
realisation of the noise. The whole time- dependent distribution of SCP-positions can 
also be obtained numerically. When the shock is localized on the left (final state), the 
width of this distribution stabilizes to a small value. The fluctuations of the shock 
velocities in different realisations of the transients states imply that the merging event 
does not occur at the same time. Hence quite large variations in the SCP positions from 
one realization to another one. Nevertheless, as we shall see, the average SCP motion 
can be predicted with high accuracy. 

6.3.1. Mean velocity of a pair of domain walls The velocity of a pair of walls can 
be computed by considering the pair as a whole (see sketch on figure H]) neglecting its 
internal structure. On each side of the "black box" surrounding the pair of shocks, 
densities are supposed to be adapted. Factorized states are again assumed on each side 
of the walls. 

We refer to quantities on the left (resp. right) of the walls with a superscript '-' 
(resp. '+'). If we note 

jf(x)^jl(a;)+ji(x) (20) 

the total horizontal flux in the left and right domains, and 

pf{x)^pl{x)+pi{x) (21) 

the total densities, then the wall velocity is given 

Atot Atot 

„tot „tot ^ ' 

p+ P- 
due to mass conservation. In case of the left pair of walls, the densities to the left of 
the walls are a^ = 1/3 and a^ = 0.2, and the densities on the right of the walls are 
pI = 0.598 and pi = 0.426. This yields a velocity 

plil - pI) + piil - pi) - a^n - a^) - a^n - a^) 

V = — ^*^ ^*^ ^*[ ^ ^ = 0.209 (23) 

pi + pi ~ a^ — a^ 

In the same way, the velocity of the right pair of walls is ^ = —0.260. This is in 
excellent agreement with the measured velocities of the second class particles, 0.209 and 
—0.261 respectively (from figure [3]-e). 

After the merging of both pairs of walls, the resulting pair of walls separates a 
right domain with densities 1-/3^ and 1 — P'^ respectively on the upper and lower lane, 
and a left domain with densities a^ and a^ respectively on the upper and lower lane. 
The predicted drift velocity of the pair of shocks is now —0.02745, and the shocks move 
to the left until they reach the boundary. Again, this velocity can be measured on the 
simulations (see fig. [3]-f), and we find — 0.02745 ±10~^, in agreement with the theoretical 
prediction. 
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Figure 4. Sketchy representation of the domain wah picture for two strongly coupled 
lanes. The walls on both lanes move together, though they are not necessarily exactly 
at the same location. Ignoring the details of the structure in the vicinity of the shocks, 
one can draw a black box around it. On each side, the densities on both lanes are 
supposed to be adapted. 



6.3.2. Structure of the shock region In the previous sections, we considered the pair 
of walls as a whole around which we can draw a black box. In this section, we would 
like to calculate, in the spirit of domain wall calculations, the internal structure of the 
shock region, and in particular the distance between the shocks. 

Outside from the pair of shocks, density profiles are flat, with adapted densities, so 
that the net vertical flux is vanishing. By contrast, in the region between the shocks, 
density profiles are x-dependent. The densities on both lanes are not adapted, and as a 
result there is a net flux from one lane to the other. 

We refer to this region as the " shock region" in the remaining of the paper. Due 
to their intra-lane dynamics, shocks would like to separate but the vertical flux acts as 
glue to keep them together. First we consider the case V = 0, i.e. the shock pair does 
not undergo any drift. 

Figure |5] gives a sketch of the simplest idealized structure that we expect. In a 
stationary state, the sum of the horizontal fluxes must be conserved through the left 
discontinuity. As j^ and p^ are continuous in x^, j^ must also be continuous in x^. 

It should be noted that, under the coarse grained description that was used in the 
previous section, the horizontal flux on each lane was discontinuous through the shock 
pair. Now that we observe the shocks at a smaller scale, for which the inner structure 
of the pair is explicitly described, only the density profiles are discontinuous, while the 
horizontal flux proflles must be continuous on each lane. Due to this zooming procedure, 
we recover the same level of description as the one that was used for weak coupling in 
1201. and thus the calculations are similar in the case V = 0. 
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Figure 5. Sketchy representation of tlie structure of a pair of walls. We assume that 
in a first approximation, we can represent the walls as real discontinuities, separated 
by a varying density profile p. 



The only way to obey the continuity constraint for j^ is to have p^{x'^_ 



1-pI. 



As 1 — pi and p_ are not adapted, a vertical flux is initialised in the shock region. In the 



same way, we must have p^[ 



Xj 



1 — pi on the left of the right discontinuity. Between 



x^ and x"^, the densities p^ and p^ must obey the coupled differential equations ( ffllSll 



with the boundary conditions 



pHx-) 

p\xZ) 

pHO 

p\xl) 



1-p 

i 



,t 



t 



P^ 



(24) 
(25) 
(26) 

(27) 



Only two of these boundary conditions are necessary to calculate the solution 
numerically starting for example at x = x"^. Of course, the solution that is obtained 
depends on the coupling strength. Matching the other boundary conditions will give 
the relation between the coupling strength and the width of the shock region x^ — x'^. 

If ^ 7^ 0, the density profiles inside the shock region are not determined anymore 
by the stationary coupled equations ( 17f8|) . Still, we assume that the shock region has a 
structure which is preserved during the drift motion. This validity of this assumption 
will be discussed later. The assumption implies that all points of the structure must be 
translated with the same velocity. 

Between the shocks, the density can be expressed by 



p^(x) 
p\x) 



p^{x 
p^{x 



Vt) 
Vt) 



(28) 
(29) 
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where the new p^ and p^ functions have to be determined. The mass conservation 
equation equivalent to fl5]|6l) but including now the time dependence is 



p^{x - Vdt) - p^{x) 
p^{x - Vdt) - p^{x) 



j\x --)- j\x + -) + h\x) - h\x) 



dt 


(30) 


dt 


(31) 




(32) 



Again, introducing a new space variable y = x/L, we find for the continuous limit 
when L becomes large 



{-V + 1 - 2p^) dyp^ = u^Lp^ (1 - p^) - u^Lp^ (1 - p^) 
{-V + 1 - 2p^) dyp^ = io^Lp^ (1 - p^) - oo^Lp^ (1 - p^) 



(33) 
(34) 



The boundary conditions for these coupled differential equations are obtained from 
the relations 



V 



p^x-) [1 - 


p^ix"^] 


]-pi 


;i - pi 




-pi 




pi 


)-pi 


pt(x-) - pi 
-p^ixl)[l-p^ixl)_ 




pI 


)-pI 


pi- p^ 

+ pi 


'i-pi_ 


-pi 


1 - pi 


i - pi_ 








pl + pi 


--pi- 


-pi 







and from the fact that p'^{x^) = pi and p^{x^. 



P\ 



Again, we have 4 boundary conditions for two coupled first order differential 



equations, and the overdetermination will allow to determine the width x^ 



X 



of 



the shock region (for a given coupling strength). When the coupling rates are of the 
same order as the particle hopping - i.e. of order 1 -, then the shock region x"^ — x^ is 
finite and the dynamics of the shocks are strongly coupled. 

This is in contrast to the case of weak coupling, for which the coupling strength 
scales as 1/L, and the distance x^ — x^ is thus of the same order as the system size. In 
the latter case, already discussed in detail by Reichenbach et al [21], the density profile 
becomes non-uniform in the bulk. 



6.3.3. Distance between the shocks of a pair We have measured the mean distance 
between the shocks in the stage corresponding to fig[3]-c, i.e. we address the time interval 
where only one pair of shocks is present in the system, which has not yet reached the 
left boundary. 

As the merging event does not occur at the same time in all realizations, the 
measurement was done between ti and t2, where ti and ^2 are realization dependent, tl 
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is the time at which the left SCP reaches its rightmost position, and t2 is the time at 
which the left SPC is located at a distance of 50 sites from the left boundary. 

We measure that (ti) = 2371 ±1.4 (with a standard deviation equal to 420) and 
(^2) = 16333 ± 8.7 (standard deviation 2745). The mean distance between the SCPs 
during this time interval ti < t < ^2 is (rf) = 4.44 ± 0.015. However, this distance 
fluctuates : the standard deviation for d is 4.8. These fluctuations can be due to two 
combined effects. First, the dynamics of the SCPs are such that SCPs perform non local 
jumps, which can be quite large. Second, the distance between the shocks of the pair 
may also fluctuate in time. To say it otherwise, the variations of the distance between 
SCPs can be due to real variations of the distance between the shocks of the pair, and/or 
to a lack of precision in the tracking of the shocks. 

In order to distinguish between the two effects, several types of tracking of the 
shocks could be compared. Different dynamics for the SCP can be defined. An 
alternative shock-tracking method uses the fact that the net vertical flux is not vanishing 
between the pair of shocks. We leave this comparison for future work. 

6.4- Final stationary state 

For the set of parameters of table [1], the final stationary state is characterized by the 
localization of the final pair of shocks near the entrance boundary. From figure [3]-f, we 
get the final SCP positions, resp. 4.9 and 7.9 on average on the upper and lower lane. 

These values could be compared with the localization lengths for the decoupled 
lanes, which are equal to ^ = In ( ^k~^: j . With the parameters of table [1], we find 
0.426 for the decoupled upper lane (wall near the entrance) and 0.272 for the decoupled 
lower lane (wall near the exit). 

It means that here, the localisation length of the pair of shocks is much larger than 
what it would be for a single shock in the decoupled lane. This is not surprising, as 
the shock on the lower lane is pulling the pair towards the bulk. Note however that we 
cannot exclude that, once again, the higher value that we find for the localization length 
is an artefact of the tracking by SCP. Another way of measuring the localization length 
would be through the density profile. An exponential fit, which is in good agreement 
with the simulation results, leads to a localization length equal to 7.1 and 7.6 on the 
upper and lower lane. 

In general deviations from the exponential profiles may occur since the applied 
domain wall picture does not consider the internal structure of the shock. Surprisingly 
for the present choice of parameters the deviations from the exponential profile are 
rather small, although the width of the pair of shocks is of the same order of magnitude 
as the localisation length. 
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value 
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value 


at 


1/3 


P' 


0.2 


a^ 


0.2 


/3^ 


1/3 


u^ 


0.2 


UJ^ 


\u^ 


A 


1 

2 







Table 2. Set (II) of parameters. Boundary conditions are adapted on both sides. 

7. Delocalization of the shocks 

The set of parameters of table [T] is such that shocks are localized at the left boundary 
in the final stationary state. Some other behaviors can be observed, for example, in the 
special case where /3t = a^ and [i^ = a^, the drift velocity of the final pair of shocks is 
exactly zero. Thus it undergoes a non biased random motion forever, and the resulting 
density profile is linear, as can be seen on fig. |6] for the parameters of table |2l More 
precisely, the linear density profile is exactly a linear interpolation between the densities 
imposed by the reservoirs, as this can be predicted by the domain wall picture. 

8. Transition from weak to strong coupling 

For the set (I) of parameters that we chose, the shocks in the decoupled system are 
located at opposite ends in the system, while coupling tends to keep the two shocks 
together. If the coupling strength is increased from zero, the first effect dominates at 
first and shocks remain at each end of the system, until the "glue" between the shocks 
is strong enough to keep them together. Figure [7] illustrates such transition from weak 
to strong coupling. In the upper lane, the shock is always localized on the left of the 
system. In the lower lane, there is a transition from a right to a left localization. 



9. Localization and adaptation 

In the whole paper, we have considered only adapted boundary conditions. In this 
section we briefiy discuss the physics of the model for more general model parameters. 
Non adapted boundary conditions create adaptation layers which interpolate between 
the non adapted densities imposed by the reservoirs, and the bulk adapted densities. 

In the case of weak coupling (exchange rates scaling as 1/L), adaptation layers 
extend over a non vanishing fraction of the system. By creating such layers on both 
sides of the system, it is possible to localize shocks in the bulk. This phenomenon is 
similar to the localization observed in one-lane ASEP models without mass conservation 
[m [131 lini [m [12]. Another way to characterize this regime would be to say that the 
width of the shocks scales as the system size, and that the density variations that we 
observe are actually the internal structure of a bulk localized shock. It is not possible 
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Figure 6. Average density profiles at different times, for tire parameters of table [2] 
and for a system size L — 100. Dashed (resp. solid) lines refer to the upper (resp. 
lower) lane. Fig (a) : density profiles at time (most upper and lower curves, black 
online) and 20 (red online). Fig (b) : density profiles at time 20 (black, red online) 
and time 200 (grey, green online). Fig (c) : density profiles at time 200 (grey, green 
online) and time 4000 (black, blue online) . Thin grey lines give density profiles at some 
intermediate times. The final stationary profile is linear and interpolates between the 
densities imposed by the reservoirs. 
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Figure 7. Final stationary density profiles for the parameters of set (I), except for 
the parameter lo^ which is varied from to 10. Here the system size is L = 100. Solid 
(dashed) lines correspond to the upper (lower) lane. 



anymore to draw a black box around the shocks, and no shocks dynamics as the one 
described in this paper can take place. 
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Figure 8. Schematic illustration of the different regions of the system in the case of 
strong coupling. Grey (green online) areas represent the adaptation regions. The pair 
of shocks (represented as dots in the bulk) is surrounded by a dashed box. 



In the case of strong coupling, the different capacities of the boundary reservoirs 
lead to an adaptation layer of a given size lau which is determined by the amplitude of 
the exchange rates, and which is independent of the system size - at least if L ;:^ lai. The 
different regions of the system are illustrated in Fig. [81 If the pair of shock is located 
in the bulk of the system, a black box can be drawn around it, and the shock velocity 
can be determined by ( 122|) since the densities on both sides of the shocks are adapted 
(and thus uniform in space). No localization is possible in this bulk region. Near the 
boundaries however, highly non-trivial dynamics is observed. If the shocks enter the 
shaded adaptation layer, the velocity of the shocks become position dependent. This 
implies e.g. that the shock might be unable to reach the boundary sites of the system, 
instead it would be reflected somewhere inside the adaptation layer. Then in the final 
state, the pair of shocks would not be localized exactly on the boundary, but at a finite 
distance from the boundary. Thus localization would also take place, but not in the 
bulk. Besides, the width of the pair x'^ — x"^ would not necessarily be small in front of 
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the adaptation layer width, so the inner structure of the pair would also be modified 
in the adaptation layer, and the shocks on each lane could be localized at different 
positions. We are not able at this stage to predict quantitatively the final localized 
state, but we expect the overall picture to be correct. 

Now we could wonder what would happen for exchange rates that scale as 1/L" 
{a < 1). We expect in this case the adaptation layers to have a width increasing with 
L. However, they should occupy a vanishing fraction of the system when L becomes 
large. Thus, some shock dynamics would be possible in the bulk, and localization would 
occur only in the adaptation layers on each side of the system. To say it otherwise, the 
extension of the shock pair could become very large, but not enough to prevent some 
motion in the bulk. The transition from weak to strong coupling can be investigated 
more systematically by applying a pertubative approach in 1 — a. 

10. Summary and conclusion 

Stochastic models of self-driven particles describe the behaviour of many transport 
processes on different scales. The self-driven particles often use one-dimensional 
structures as tracks, e.g. biopolymers or, at a larger scale, highway lanes. Although 
being one- dimensional, often several tracks are offered to the dynamic particles. 

In this work, we consider the case of two-lane systems with particle exchanges 
between the lanes. We focus on the case of strong coupling, i.e. particle exchange rates 
are of the same order as hopping rates. In contrast to previous work we concentrate on 
the transient regime rather than the stationary behaviour. We show that, in the strong 
coupling limit, the transient regime can be fully understood in terms of a domain wall 
picture, i.e. in terms of the dynamics of pairs of shocks. For concreteness we mainly 
discussed two sets of parameters which are chosen in order to illustrate the generic 
behaviour of the model in the case of strong coupling. 

Our results show that the transient regime can be divided into two main stages. 
The first stage is a very rapid local adaptation between the lanes. This stage occurs on 
a time scale which is determined by the amplitude of the exchange rates, while the bulk 
densities at the end of this stage are determined by the ratio of the exchange rates. After 
this first rapid stage, the exchange fiux between the lanes is zero almost everywhere in 
the system, i.e. everywhere except in some layers near the boundaries, as this will be 
discussed below. 

In the second stage, we have shown that the dynamics is dominated by the dynamics 
of several pairs of shocks, which undergo drift motion, merge when they meet another 
pair, and drift again until the unique final pair reaches the boundary where it stays 
localized. In this slow relaxation stage, the pairs of shocks have a dynamics very similar 
to the one observed in mass conserving systems. The non-conservation of mass in each 
lane is concentrated between the shocks of a same pair, i.e. the vertical fiux is zero 
everywhere in the bulk except between the shocks of a same pair, where it acts as a 
"glue" that keeps the two shocks together. In this stage, in order to track the shock 
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positions, we have generalised the dynamics of the so-called second class particles to 
two lane systems. These results underline the flexibility of the domain wall approach 
for one- dimensional transport. 

Near the boundaries, the densities imposed by the reservoirs can prevent local 
adaptation, at least in finite size layers called adaptation layers. In this paper, we have 
chosen the parameters of the model such that the boundary reservoirs are adapted, i.e. 
the net-fiux between the lanes is zero also near the system boundaries, and there are 
no adaptation layers. This choice allows us to characterise the relaxation of the system 
purely via the domain wall dynamics. However, our results could be easily extended 
to non adapted reservoirs as, in the case of strong coupling, adaptation layers have a 
constant width - and thus represent a vanishing fraction of the system size when the 
system size increases. Then the dynamics of the shock pairs is modified only in the 
adaptation layers, but remains unchanged in the bulk. 

In the weak coupling limit, adaptation layers invade the whole system, and bulk 
localization of the shocks can be observed [201 EI] • Here, for strong coupling, adaptation 
layers can prevent the shock pairs from reaching the boundaries, and can localize the 
pairs at a fixed distance from the boundary - but no localization is possible in the bulk. 

While our focus was here on the domain wall dynamics in the transient regime, [33] 
have developped a systematic approach to determine the phase diagram for the final 
stationary state, in the case of adapted boundary conditions, which should shed more 
light on the final coupling state between the walls. 

In principle, the domain wall approach should allow to determine the stationary 
phase diagram in the strong coupling case, though we chose to address only the non- 
stationary dynamics here. This is a supplementary illustration of the predictive power 
of the domain wall picture, even in relatively complex systems. 
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